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1 Introduction 

The problem of distinguishing on the quantum level between integrability and chaos 
of classical systems is a recurrent, if not the principal, topic of the quantum chaos 
theory. The usual definitions or signatures of the chaotic motion based on such 
phase-space notions like Lyapunov exponents or various degrees of ergodicity lack 
sense on the quantum level due to the absence of the proper notion of the phase 
space in quantum mechanics. Instead, in search of criteria of quantum chaos one 
should resort to purely quantum characteristics of a system, like e.g. its spectral 
properties. 

A quarter of a century ago, Bohigas and coworkers PQ proposed a character- 
ization of quantum chaos based on the statistical theory of spectra. According 
to their hypothesis most quantum systems whose classical limit is chaotic display 
universal spectral fluctuations determined by the Random Matrix Theory (RMT) 
[2]. Generic classically nonintegrable systems exhibit thus level repulsion i.e. the 
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probability of finding two adjacent energy levels tends to zero with the difference 
between their energies. 

A vast numerical and experimental evidence [3, 4] in favor of this hypothesis 
was collected during last twenty years. At least two general strategies providing 
theoretical arguments supporting the Bohigas-Giannoni-Schmidt conjecture were 
developed. The first one employs statistical mechanics of a fictitious gas of eigen- 
values undergoing parametric dynamics, where the role of time is played by a control 
parameter controlling the transition from an integrable to a chaotic system [3l H] . 
The other approach takes its roots in the semiclassical quantization via classical 
periodic orbits pioneered by Gutzwiller [2] and extends ideas from the theory of 
disordered systems to dynamical systems [6]. In both approach, at some stage one 
invokes some statistical hypothesis which justifies ascribing statistically inferred 
properties to an individual system which is the object of actual numerical or ex- 
perimental examinations. On the other hand the proofs of nonintegrability on the 
classical level are usually based solely on numerical investigations. We do not know 
of any example of a model system for which there exist analytical proofs of the 
classical nonintegrability on one side and of the repulsion between quantum levels 
on the other. The aim of the presented investigations is to provide such a model. 
In the present paper we will concentrate on the classical side of the problem and 
show its classical nonintegrability. 

In a series of papers EJ [U] one of us proposed a class of models taking their 
origins in atomic physics and quantum optics in which dynamical variables were 
elements of a compact semisimple Lie algebra in some particular irreducible rep- 
resentation. The classical limit was attained by going with the dimension of the 
representation to infinity. For Lie algebras with the rank larger that one (e.g. for 
SU3 algebra) there are more than one, 'natural' ways of performing this limiting 
procedure. In the effect there are several inequivalent classical corresponding clas- 
sical systems, differing e.g. by the dimensionality of the classical phase space. The 
interplay between the number of degrees of freedom and the dimensionality of the 
space is crucial for the (non-)integrability of a classical Hamiltonian system. A 
particular quantum system can be classically integrable or not, depending on the 
way the classical limit is approached. The story becomes interesting if, basing 
on the above observation, one can say something about purely quantum features 
(e.g. spectral properties) of the quantum system in question. The affirmative reply 
based on numerical investigation of spectra and classical characteristics of chaos, 
was given in [8]. 

In the present paper we examine again the above models in order to prove in 
an analytical, rather than numerical manner the nonintegrability. Admittedly it 
is a rather minimalistic goal. Nobody claims that mere nonintegrability of the 
classical system is sufficient for the repulsion of the quantum energy levels, we 
believe that the Bohigas-Giannoni-Schmidt conjecture gives correct predictions only 
for 'sufficiently chaotic' systems. Nevertheless the nonintegrability is a necessary 
prerequisite and we decide to pursue such a modest goal of showing it for the 
considered system. 
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2 Nonintegrability of Hamiltonian systems 



In this section we present a brief introduction to the Morales-Ramis theory of 
non-integrability of Hamiltonian systems \10\ [TT] . The aim is to give some basic 
definitions and theorems (without proofs) and to present general scheme which 
should be followed when one wants to prove non-integrability. The whole concept 
is based on the differential Galois theory (here good references are [121 HU [H] ) , so 
we also give some key ideas of it. 

2.1 Differential Galois theory 

Differential Galois theory (a.k.a. Picard-Vessiot theory) is an analogue of the clas- 
sical Galois theory (which deals with algebraic equations), for linear differential 
equations. The main object of the theory is here a differential field K i.e. an 
algebraic field equipped with a differentiation, i.e. a linear mapping ' : K — > K 
satisfying the Leibniz rule, (fg)' = fg + fg' ■ 

The object of our interest is a linear differential equation: 

+ dn^-V + ... + ai y' + a y = 0, (1) 

where the coefficients a% are elements of K (it is good to think K = C(x), the field 
of rational functions). The natural question which arises in the differential Galois 
theory is: can one solve the equation (TjQ) using some "elementary" operations? Of 
course usually the solution of (pQ) is not contained in the field K and we have to 
extend if to a larger one. The smallest extension of K which contains all solutions 
of ([I]) is called the Picard-Vessiot extension. It exists always when the field of 
constants of K is algebraically closed |12j . 

To be more precise we have to define what we mean by "elementary" operations. 

Definition 1 Let K be a differential field. We say that the equation (Tip is solvable 
in the Liouville functions category when the Picard- Vessiot extension of K can be 
obtained from K in a finite number of steps each one being an extension of K by 
adding: 

1. on which is algebraic over K, 

2. on which is such that cJ i G K , 

3. on which is such that € K. 

The steps 1., 2., and 3. correspond to the operations of taking roots of polynomial 
equations, integration, and taking the exponent of an integral — the natural op- 
erations when one solves a differential equation. The last important definition is 
that of the differential Galois group: 

Definition 2 Let L D K be a Picard-Vessiot extension of K for the equation (TT]). 
The differential Galois group Gal(L D K) of the extension L D K is the group 
of all differential automorphisms ( algebraic automorphisms that commute with the 
differentiation) of L which are identity on K. 

Let y\ be a solution of (pQ) and a £ Gal(L D K). Then: 

aiy^ + a^y^ + ... + a iy [ + a yi) = a(0) = (2) 
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Using properties of a, namely the fact that it is a differential automorphism, we 
obtain: 



{yi) (n) + a n - 1( x{yx)( n - X) + . . . + axa( yi )' + a a{ yi ) = 



(3) 



The last equality says that every element of Gal(L D K) gives a solution of (TjQ) 
when acting on a solution of ([T]). Let {yi,y2> • • • ,Vn} be the fundamental set of 
solutions of dU). This means that any solution of ([T]) has the form y = Yli=i a iVi 
where 04 are some constants, and 



From the last equality we see that any a £ Gal(L D K) is completely determined 
by its action on solutions of (pQ), i.e. 



where a, are constants. As a result we can represent the differential Galois group 
as a subgroup of GL(n, Ck) where Ck - the field of constants of K. In fact this is 
an algebraic subgroup of GL(n, Ck) [IS], so in particular a Lie group. Now we can 
formulate the main theorem [T2] 

Theorem 1 The differential equation |7]j is solvable in the Liouville category if 
and only if the corresponding differential Galois group is solvable. 

Thanks to the fact that Gal{L D K) is a Lie group it is enough to check whether 
the Lie algebra of Gal{L D K) is solvable. 

2.2 Morales- Ramis theory - the general scheme 

The Morales-Ramis theory is a powerful tool for checking (non-)integrability of 
Hamiltonian systems. By an integrable Hamiltonian we understand here one which 
admits enough number of functionally independent, involutive with respect to the 
Poisson bracket integrals of motion (this number should be equal to the number of 
degrees of freedom of our system). Establishing the nonintegrability of a Hamil- 
tonian system is thus equivalent to proving the nonexistence of the appropriate 
number of integrals of motion. 

Let us shortly outline the logic of such proofs. For a general system of nonlinear 
differential equations a direct count of the number of integrals of motion is difficult 
- in principle there are no known methods of achieving the goal. It is, however, 
clear that integrability of a nonlinear system should be, in some way, inherited by 
its linearized version. Now, for linear equations we have powerful methods based 
on the differential Galois theory outlined above which can be used to relate the 
(non)-integrability to properties of the differential Galois group. Reversing now the 
argument we see thus that if by analyzing the differential Galois group we are able 
to establish nonintegrability of the linearization, we will prove the nonintegrability 
of the full, nonlinear system of equations. An additional bonus is provided by the 
Hamiltonian character of the system which simplifies the structure of the Galois 
group of its linearization. 





n 




(5) 



i=l 



4 



To make the above outlined idea precise and workable we need some concepts 
and facts. Let (M, uj) be a symplectic manifold, i.e. a; is a closed (du = 0) and non- 
degenerate two-form on M. The Hamilton equations corresponding to a Hamilton 
function H have the form 

x = X H (x), (6) 

with Xh defined via ix H u := u(Xh, •) = dH. A linearization of © is now achieved 
by defining the variational equation (VE) along a particular non-equilibrium in- 
tegral curve r : x = (j){t) of ([6]). By considering a solution of © in the form 
x' = <p(t) + x(t) an d retaining only the linear terms we obtain the familiar result: 

X = X' H (cP(t)) X . (7) 

For further considerations concerning additional possible simplifications and reduc- 
tions of (J7|) it is worth treating the derivation of it on a slightly more formal level. 
We first observe that we can restrict the tangent bundle TM to T obtaining the 
vector bundle TM\-p over T. We define operator D on TM\y to be the Lie deriva- 
tive Cx H restricted to TM\-p. More precisely to compute &x H Y we extend Y to Y 
on a neighborhood of T, compute Cx H Y, and restrict the result to T. Operator D 
inherits the properties of the Lie derivative, in particular D{fY) = f'Y + fD(Y). 
The variational equation along V is simply DY = 0. It can be easily checked that 
choosing a suitable basis in TM |r it can be written in the form (|7|). 

Now we want to explore the above mentioned idea that if the system © is 
integrable, then the system ([7]) is integrable as well. To this end we use the Ziglin 
lemma stating that with every first integral / of the system ([6]) we can asso- 
ciate a first integral f° of (J7|). Moreover if fx, fe, . . . , are involutive, functionally 
independent first integrals of © then the corresponding functions f°, ...,/? are 
involutive, functionally independent first integrals of (|7|). The Morales- Ramis' idea 
was to investigate which restrictions on the differential Galois group of the vari- 
ational equation are imposed by the complete integrability of the system ([6]). It 
turns out PU EE] that: 

Theorem 2 Assume that a complex analytic Hamiltonian system is integrable in 
the meromorphic function category, then the identity component of the differential 
Galois group of the corresponding variational equation is abelian. 

The complete integrability indeed imposes a very strong condition on the differen- 
tial Galois group of VE. To use this theorem effectively we notice that if we know 
k involutive integrals of motion of the system © we can reduce the dimensionality 
of the system - this is a familiar procedure known from standard classical me- 
chanics. More precisely the involutive integrals of motion determine k commuting 
Hamiltonian vector fields which, in turn, define an isotropic0 subbundle F of TM. 
We can now perform a symplectic reduction. We have a well defined symplectic 
form on F u /F = F N and, further, we can restrict the operator D to F because 
DY G F^-^ if Y € F^ . After the reduction we obtain from (|7|) a 2(n — k) - dimen- 
sional system called the normal variational equation (NVE) [10|, [TT] . It is proved 
in [10] that if the differential Galois group of VE is abelian then this also true for 
the differential Galois group of NVE. Moreover if there is only one missing first 

1 A subspace W of a symplectic space V is isotropic if and only if W C W i.e. W is a subspace of 
its orthogonal complement in the sense of the symplectic form to. 
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integral of motion then VE can be reduced to a 2-dimensional linear matrix dif- 
ferential equation equivalent to some second order linear differential one. Now the 
investigations of the Galois group can be performed in an algorithmic way. From 
the general theory (see Section 12. ip we know that the Galois group is a subgroup of 
GL(2, C). In fact we can restrict our search to subgroups of SL(2,C) [12]. Indeed, 
by a change of the dependent variable z(t) = exp (^A(tf\ y(t), with A'(t) = a\(t) 
we can eliminate the first-derivative term from the equation y" + a\y' + a^y = 
obtaining z" + b§z = with bo = oo — \oq — \a' Q without spoiling such properties of 
the coefficients like meromorphicity or rationality. It is a matter of a short calcu- 
lation to show that the Wronskian, W = z\z' 2 — z' 1 Z2, of two solutions Zi, Z2 of the 
new equation is a constant function, a non-zero one if z\ and zi are independent. 
Hence, from the definition of the differential Galois group we have a(W) = W 
for its arbitrary element a. On the other hand, by a straightforward calculation, 
a(W) = det(a)W, hence det(cr) = 1. 

Among four possibilities allowed [12], i.e. the differential Galois group being 

1. a finite group: the tetrahedral group, the octahedral group or the icosahedral 
group, 

2. the group of matrices conjugated to the subgroup 





' c 
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c- 1 


7 


c- 1 


i.e. matrices of the form AX A' 


-1 

1 


where A 



where A is a fixed element of SL(2, C) and 



X varies over the subgroup, 



the group of triangulizable matrices, i.e. matrices conjugate to the subgroup 

c d 



of triangular matrices, 



c 



i 



and 



4. the whole ST(2,C), 

only in the last case it is not solvable. A practical way of establishing the relevant 
case for a particular equation is provided by the Kovacic algorithm [16] which can 
be used to determine the differential Galois group upon analyzing poles of the 
coefficient bo. Of course, if the group is not solvable then it is not abelian either, 
hence our system is not integrable. 



3 Classical and quantum chaos of 5U3 systems 

In [8] a special class of systems having a compact phase space on the classical level 
and, consequently, a finite-dimensional Hilbert space in the quantum setting, was 
investigated. The classical limit is approached by increasing the dimension of the 
Hilbert space. 

As an example let us consider a collection of N atoms interacting resonantly 
with the electromagnetic radiation. Usually, due to imposed resonance conditions, 
only a finite number n of energy levels of each atom is involved in the interaction. 
Transitions from the level \l) to the level \k) of a single atom are described by the 
operators sjy = \l){k\, which span the defining representation of the Lie algebra 
gl n (C) in the n-dimensional Hilbert space spanned by the states of an atom. In a 
system of N atoms confined to a small volume in which they feel the same field 
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amplitude the transitions are described by the operators Ski = Yl a =i s kv wnere s ki 
acts as Ski on the levels of the a-th atom and as the identity on the rest. Clearly Ski 
span the N-th tensor power of the defining representation of gl n (C) relevant for a 
single atom. The resulting representation is clearly reducible. Various preparation 
of the initial state of the whole system of atoms determine its relevant irreducible 
components. 

Typically the number of atoms is conserved, so N = Y2i=i ^ii 1S a constant of 
motion, and we can restrict the considerations to sl n (C). The observables of the 
considered model are constructed as polynomials in the generators. Since they have 
to be hermitian we finally focus our attention on the g = su n algebra and G = SU n 
Lie group. 

The dynamics of the observables is governed by Heisenberg equations of motion 
generated by the Hamilton operator of the considered system. The classical limit 
becomes relevant when we increase the number of atoms and ask questions about 
such quantities like e.g. energy or polarization per one atom. Formally it consists 
of putting N —* oo. We expect that in the limit the generators Ski, after appropri- 
ate scaling (e.g. by the number of atoms), are mapped into classical functions on 
appropriate phase space in such a way that the Heisenberg equations of motion are 
mapped to classical Hamilton equations. In this way the 'Dirac quantization' pro- 
cedure requesting correspondence between commutators of observables and Poisson 
brackets of the corresponding phase-space functions ('classical observables') is ob- 
served. 

Increasing the number of atoms N results in the growing dimension of the 
largest irreducible component of the constructed representation. The construction 
of the classical phase space is achieved by the following limiting procedure. For a 
irreducible representation of G in a vector space V we consider its projective variant 
i.e. the action of the group on the projective space P(V) given by g ■ [v] = [g ■ v] for 
g £ G, v £ V and [v] £ P(V) - the ray through v. It is known [IT] that the orbit of 
G through the point [v] € P(V) corresponding to the highest weight vector, i.e. the 
common eigenvector of all Sa annihilated by all Sy with i < j, is endowed with 
a natural symplectic structure. The projective orbits through the highest weight 
vectors can be mapped on the orbits of the coadjoint representation of G [IB] i.e. the 
representation of G on the dual space 0* of the Lie algebra g with the corresponding 
symplectic structure known as the Kirillov-Kostant-Souriau form. Coadjoint orbits 
are thus good candidates for classical phase-spaces. Each irreducible representation 
and each above defined orbit is uniquely determined by the highest weight vector 
V or, equivalently, by the corresponding eigenvalues of Sa, the number of which 
equals to the rank of the group. Alternatively and equivalently, to identify an 
irreducible representation we may use independent Casimir invariants - elements 
of the enveloping algebra of g. From the definition they commute with all elements 
of g, hence for each irreducible representations they are constant multiples of the 
identity operators. The values of these constants identify a representation and 
consequently also an orbit. 

The symplectic structure on coadjoint orbits can be also obtained from a nat- 
ural Poisson structure on the linear space g*, so called Lie- Poisson structure (see 
next Section). Its symplectic leaves, i.e. manifolds on which the Poisson bracket 
determines a non-degenerate two- form, are exactly the orbits of the coadjoint rep- 
resentation. Finally thus we can identify the classical phase space with a symplectic 
leaf of the Lie-Poisson structure on g*. 
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For each element of the sequence of irreducible representations with growing 
dimensions we obtain thus a unique well defined classical phase space. After an 
appropriate scaling by the volume of the orbit [Hi] we obtain in the limit the desired 
phase space of the limiting classical system, dynamics of which is connected with the 
quantum one via Dirac's correspondence. The general idea of this construction goes 
back to Simon |20j : see (9j[T9] for the setting relevant for the present considerations. 

The above presented construction of the classical limit is purely geometrical. 
To make it more appealing form the physical point of view let us observe that the 
final result can be also obtained by treating expectation values of the quantum 
observables as classical phase-space functions in the limit of vanishing Planck con- 
stant [8] . The relevant expectation values are calculated for appropriate coherent 
states of the group G - these are the states 'most classical' from the point of view 
of uncertainty principle, hence the best approximations to the classical description 
of a system [2T] . 

The symplectic leaves of the Lie-Poisson structure may have different dimen- 
sions. In the construction outlined above the dimensionality of the resulting phase 
space may thus depend on the chosen way through the sequence of irreducible rep- 
resentations. In the simplest non-trivial case of three-level atoms the algebra of 
observables is spanned by eight quantities, the generators of the Lie group SU 3 . Ir- 
reducible representations of SU 3 are indexed by two independent quantum numbers 
- the weight of the highest- weight vector [22\ - which determine also the dimension 
of the representation. In effect there exist two inequivalent ways to the classical 
limit resulting in a six- (in a generic case) or four- (in a degenerate case) dimen- 
sional classical phase space. This observation was a basic point of the paper [8], 
where it was shown that the dimensionality of the classical space determines not 
only integrability properties of a specific class of classical Hamilton functions, but 
also some statistical properties of spectra of the quantum Hamiltonians for which 
the classical system in question is the classical limit outlined above. Both investi- 
gations of spectra on the quantum level and integrability properties on a classical 
one were performed numerically. We are now in position to prove analytically the 
non-integrability for a concrete member of the considered class. 

The Hamiltonian we consider is quadratic in the generators Sij, 

H = 3(S? 2 + S 2 21 ) + 15(S 13 S 32 + S 23 S 31 ) (8) 

(see remarks in [71 [8] for the possibilities of experimental realizations). It is easy 
to show that [H, Y] = 0, where Y = Su + S 22 — 2S 33 . The commutation relation 
survives the classical limit providing thus an integral of motion. As a result we 
obtain a classical Hamiltonian system with the Hamilton function 

H = 3(s? 2 + s li) + 15(si3S32 + S23S31), (9) 
admitting an integral of motion 

^ = S11+S22-2S 3 3, (10) 

where are coordinates on SU3 dual to . 

In the next section we explain how to obtain from Q the corresponding Hamil- 
ton equations of motion on SU3, in particular we describe the above mentioned 
Lie-Poisson structure. Since the classical phase space is either four or six dimen- 
sional the classical system ([9]) is integrable in the former and possibly non-integrable 
(if there are no other unknown integrals of motion) in the latter case. 
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4 Lie-Poisson structure on sut 



As stated in the Introduction and explained in the preceding section, we are in- 
terested in a class of systems which are obtained as the classical limits of some 
quantum systems with SU(3) symmetry [7J EJ [9], and the classical limit of such 
a system can be considered as a Hamiltonian system on SU3. It is instructive to 
consider the problem in a slightly more general setting where SU3 is substituted by 
an arbitrary Lie algebra. 

Let (g, [•,•]) be an arbitrary Lie algebra. We can equip its dual space g* with 
the canonical Lie-Poisson structure given by the Poisson bracket, 

{f,g}(x) = (x,[(df) x ,(dg) x ]), (11) 

where f,g £ C°°(q*) and (,) denotes the pairing between g* and g. The bracket 
(fTT]) is well defined because (df) x : g* — > R is an element of g** and in the case of 
finite dimensional vector spaces we have g** = g. The bracket defined in this way is 
of course bilinear and antisymmetric. It is also a differentiation i.e. it satisfies the 
Leibniz rule. It is convenient to describe the Lie-Poisson structure on g* in terms 
of a Poisson bivector, 

^ =c Nfe A ^? (12) 

where c& are the structure constants of g corresponding to some basis ei,...,e n , 
i.e. [ei,ej] = c^e^, and coordinates the vectors ei,...,e n considered 

as linear functions on g*. To show that the bivector induced by Poisson bracket 
(jlip is the same as the one given by Eq. (|12p . it is enough to check it on the 
linear functions. After completing this easy taks we may thus write the Hamilton 
equations of motion for na arbitrary function / on g* , 

& = {HJ} = r,(HJ) (13) 

The Poisson bivector r] is degenerate on g* due to the existence of Casimir 
functions which have a vanishing Poisson bracket with any function. Thus we do 
not obtain directly any symplectic structure on g*. On the other hand when we 
restrict rj to its symplectic leaves determined by constant values of independent 
Casimir functions we end up with well defined symplectic manifolds. Indeed, rj 
defines a morphism 

v i : Tq ** -» Tg*, df\ x ^X f \ x , (14) 

where Xj = {/,•}. It generates a distribution V = [j x D x , x G M, and D x = 
image(r) x ). This distribution is involutive, (i.e. [X, Y] G T>, for X, Y G T>), hence 
from the Frobenius theorem T> is tangent to some generalized foliation T . The 
restriction of rj to leaves, J- x , x G M, of the foliation J 7 is a well defined, non- 
degenerate Poisson bivector, so it defines a symplectic structure on F x . The sym- 
plectic leaves J- x are exactly the coadjoint orbits of G, and the corresponding sym- 
plectic form providing a symplectic structure is the announced Kirillov-Kostant- 
Souriau one [18] . Its explicit form can be easily deduced from the definition of rj, 
It is, however, often easier to work with the Poisson structure (lllh on the whole g 
than with its restriction to leaves and treat the Casimir functions as constants mo- 
tion determining by their initial values a manifold to which the motion is restricted. 
This is the way we will follow in our case. 
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We may now specify the above general considerations to the SU3 case of the 
present interest. To this end we have to chose some basis SU3 and find explicitly 
the bivector 77. The basis of our choice consists of the standard Gell-Mann matrices 
(see e.g. [22J; we use them in slightly different order) multiplied by the imaginary 
unit i, 



ei 



C-l 



?7 



10 
-10 

1 

-10 

% 
i 
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i 


" 




i 
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, e 2 = 


i 








, e 3 = 
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, e 5 = 











, e 6 = 








1 
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-1 






(15) 
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i 







-2t 



In this basis we have 77 
given as 



a 



Vij-g%- A g|- with the coefficients rjij 



V12 


= 2x 3 , 


V25 


= -X6, 


V46 = 


-XI, 


V13 


= -2x 2 , 


me 


= x 5 , 


V47 = 


X2, 


riu 


= ~xe, 


r]27 


= -x 4 , 


V48 = 


-y/3x 5 , 


V15 


= -X7, 




= x 5 , 


r}m = 


-X2, 


Vie 


= x 4 , 


%5 


= -x 4 , 


V57 = 


-Xl, 


Vn 


= x 5 , 


??36 


= -X7, 


V58 = 


\/3»4, 


V23 


= 2xi, 


%7 


= Xq, 


V67 = 


(-x 3 + \/3x 8 ) 


mi 


= x 7 , 


^45 


= (x 3 + \/3x 8 ), 


V68 = 


-V3x 7 , 



V78 



\/3x 6 . 



The bivector 77 is degenerate and there are two functionally independent Casimir 
functions c%, C2, (rf{dc\) = = rf(dc2)), given by c± = atr(X 2 ) and C2 = /3tr(X 3 ), 
where X is a generic matrix belonging to SU3 algebra, 



X 



ix 3 + ^§ 

-xi + ix 2 
-x 4 + ix 5 



xi + ix 2 
-ix 3 + i% 

-X 6 + X7 



x 4 + ix 5 \ 

Xq + tX 7 



Zl V3 



J 



and a, (5 are arbitrary constants. To keep the consistency with [S] where hermitian 
rather than antihermitian matrices were used to represent SU3 algebra, we choose 
a = —1, j3 = i. 



5 The non-integrability proof 

As we noticed in Section [H the SU3 Lie-Poisson structure has two Casimir functions 
ci and C2, hence the dimension of a generic leaf is six, but there are cases for 
which it reduces to four [7j. It is interesting to check whether a Hamiltonian 
system defined on the whole SU3 which possess an additional first integral besides 
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the Hamilton function itself, 'placed' on leaves of different dimension is integrable 
or not. By 'placed' we mean that initial conditions determine the leaf on which 
the time evolution takes place. Of course such a Hamiltonian system is integrable 
on four dimensional leaves (there are two first integra Is in involution), but in the 
case of six dimensional leaves an additional first integral needed for integrability 
may lack. Our aim is to prove that for some particular polynomial Hamiltonian 
systems an additional first integral is indeed missing, with ([9]) treated as a concrete 
example. 

Hamilton functions H on SU3 we are interested in are given as second order 
polynomials in the Xj coordinates given in Section 01 The particular example given 
by © takes in the new variables the forrrH, 

H = 6(xf — x\ — 5x4X6 — 5x5X7). (16) 

The resulting Hamilton equations are given by: 

^=V 9 (dH). (17) 
In the coordinates x% they form a set of eight differential equations: 

^ = 24x 2 x 3 + 30xj + 30x i ~ 30x i ~ 30x 7> 
dx\ 

~dt 

dx 3 

dt 

dX4 

~~dt 

dx 5 

dt 
dxQ 

~~dt 
dxi 

~~dt 

dxg 

~dt 



24xix 3 , 

— 48xiX2 + 6OX5X6 — 6OX4X7, 

— 30xiX4 + 30X2X5 — 12xiX6 — 12x2X7 + 30X7(X3 + V3xg), 
— 30X2X4 — 3OX1X5 + 12x 2 X6 — 12x1X7 + 30xg(— X3 — \/3xg), 
12x1X4 — 12x 2 X5 + 30xiX6 + 3OX2X7 + 30x5(— X3 + v^xg), 

12x 2 X4 + 12xiX5 — 30X2X6 + 30xiX7 + 30X4(X3 — \/3xg), 

0. (18) 



The last equations reflects the fact that (|10p is a constant of motion, since in the 
new coordinates we have simply Y = xg. As we announced we are working in the 
full su* space, we know thus two additional integrals of motion given by the Casimir 
functions c\ and 0%. In [8] it was shown that the for the classical limit of SU(3) 
systems obtained via the procedure outlined in Section [31 the values of the Casimir 
constants of motion can be parameterized by a single number q £ [0, 1], 

ci = 3(9-9 + 1), 

c 2 = ]-(-2q 3 + 3q 2 + 3q-2). (19) 

It was also shown in [8] that for q = and q = 1 the leaf on which the system 
evolves is four dimensional whereas for q e]0, 1[ the leaves are six dimensional. 



2 The change of variables from Sjj to Xi can be easily red off from (TJ]) via Xij = Sj 
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Fixing the values of c\ and C2 by choosing a particular q restricts also possible 
values of x%. A possible compatible choice is 

X8 = ^(l-2q). (20) 

To use Morales-Ramis theory we have to find some particular, as simple as 
possible, non-equilibrium solution. If we put 

Xl = x 2 = X 5 = X 6 = Xg = 0, X 4 = X7, (21) 

the system (fl"8j) reduces to two differential equations: 

^ = "60*i (22) 
= 3OX3X4. 

dt 



The system can be easily solved since we know two constant Casimir functions. We 
want to choose such values of them that the corresponding symplectic leaf is six- 
dimensional. For our choice (|2ip the Casimir functions simplify to c\ = 2(x3 + 2x|) 
and C2 = 0. Observe that a choice q = ^ which, as stated above corresponds to 
a six-dimensional phase space, indeed gives C2 = 
xg = [see (f20|) ]. We will prove that in the case q 
not integrable in the Liouville sense. 

The solution of (I22p for q = ^ is found to be 



[see (|19p ]. forcing in addition 
= i the whole system (fl"8j) is 



X3 = — tanh(— 15t) 



X4 



V2 



a/1 - tanh(-15t), 



(23) 



It defines a particular integral curve T of the full system. We can now find varia- 
tional equation along the obtained solution, 



dy 
dt 



A(t)y, 



(24) 



where the matrix A(t) is given by 



Ait) 



\ 






24x 3 





6OX4 








—6OX4 





24x 3 
































—6OX4 








—6OX4 





—30x4 


-12x 4 


30x4 











30x 3 


30\/3x4 


-12x 4 


—30x4 











-30x 3 








12x 4 


30x4 








-30x 3 











30x4 


12x 4 


30x4 


30x 3 











— 30\/3x4 



























/ 



The tangent space at each point x of curve T is D x where D x = image{j] x ). 
It can be described effectively using Casimir functions, namely V x G T X F X if 
and only if dci(V x ) = and dc2(V x ) = 0. If we denote by V x = {V x , . . . , V^ 8 } 
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the coefficients of V x with respect to the basis corresponding to x\, . . . , Xg, the 
conditions dc\(V x ) = and dc2(V x ) = take the form: 



XsV^ + x^ + V]) = 0, 
2,x\V x l + 3x3 + 24 (X 4 - Vl) + V3(xl - x\)V* = 



(25) 



Using these relations we notice that V x and V x are completely determined if 
we know V X ,V X , Vj, so we can reduce (1241) to a set of six differential equations 



dt 

where the matrix B{t) is given by: 
/ 



B{t) 



\ 






24x 3 


6OX4 








24x 3 














30x4 (x§ +2x1) 

2" 2" 

*^ 3 


-12x 4 


30x|(-4x§+x|) 
3 2 

Xg X 3 X^ 








-12x 4 


—30x4 








-30x 3 


12x 4 


30x4 





-30x 3 





30x 4 + 4% 


12x 4 


30(xf+x2x|+x|) 

_3 '',.2 









(26) 



-60x 4 \ 


30 ("-^2 ~\~ "-^*3 4 — \~ 4. ) 
3" 2" 






30x2(-4x2+x|) 

5 — r^n^ — 

Xg X 3 X^ 



/ 



The set of differential equations ( 126]) can be reduced to a normal variational 
equation (NVE) using the fact that X H = rf(dH) and X s = tf{dx & ) are 
solutions of (1261) . To perform the reduction we need the symplectic form at 
each point of the curve V. It can be obtained by inversion of rj-p which is 
the restriction of 77 to the curve T. It is possible since 77 is non-degenerate 
on r (r C F x for some x). Explicit calculations give the symplectic form 
along r as u r = co^dxi A dxj, where i and j belong to {1, 2, 4, 5, 6, 7} and the 



coefficients Uij read: 



U12 = 


X3 


2(xi-xi) 


Ul5 = 


X4 


2(xi-xl) 


^16 = 


X4 


2(xl-xi) 


UJ24 = 
^45 — 


X4 


2 a? ■ ^4 


2x 3 (x§-x|) 



a>27 

<^45 
^46 
^57 



2x§ 



2x 3 (x|-x|) 

2x 3 (x|-x|) 

— a;^ 
2x 3 (x|-x|) 



The key point now is to find a symplectic basis including X# and X ; 
set of six vector fields such that r/r ~ X# A + X$ A X$ + X A X . 
basis always exist [TO], and in our case it is formed by 



8 , i.e. a 
Such a 



X 



H 



X* 



30 X3X4 



<9 



+ 







X4 



X = -2x 



dx 5 
d 

3 dx\ 



dx 4 dx 
_d_ 

dxt 

d 

+ x 4 



X 



H 







d 



dx^ dx-j 



X* 



X 



60 X4 

rp 2 ™ 2 

2 Z3X4 



+ 



5 



dx 



d 
8x4 
X4 

I — - 

<9x 2 2x 3 



+ 







dx-j 











dx<s dxe 



(27) 
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We can express the equation (I2"B"]) in the basis (1271) as: 



^ = p- 1 (S(t)P-P) = C(t) X 
where x = P£, an d -P is the change of basis matrix: 
/ 



(2J 



P(f) 






3OX3X4 



\ 30x3X4 








X4 
— x 4 










-1 

6OX4 

-1 

6OX4 








2 2 
x 3~ x 4 

2x3X4 





2x3x4 



-2x 3 


X4 



— X4 





-1 



X4 
2x 3 
_ x & 
2x 3 



/ 



The derivative matrix P can be easily computed using ( 1221) . The final result 
reads as 



C(t) 



Thus the NVE is a simple 2x2 matrix differential equation: 



















2 \ 

















30 ^0-^4 
x 3 















































-30 + ^4- 

x 3 
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\o 











48x1 


/ 



dx 
dt 



48x1 
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X- 



Writing 



X 



Xi 

X2 



we find the corresponding second order differential equation, 

d 2 Xi 



dt 2 



- 576x3X1 



0. 



Making use of 



we get: 
d 2 X 



dt 2 



- - 144tanh 2 (-15t)x! = 



The substitution y = tanh(— 15t) transforms the equation to: 

255 (1 - y 2 ) 2 - 450 (l - y 2 ) y^ - U4y 2 f = 



(29) 



(30) 



(31) 



(32) 



The equation is fully prepared to be treated by the Kovacic algorithm. We 
are not going to describe it here (see the references [T6j and [TO] for details). 
The algorithm produces a solution if an equation is integrable in the Liouville 
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sense, and, what is more important for us, it determines as a byproduct the 
differential Galois group identifying it as one among those listed at the end 
of Section 12.21 For our equation (1521 the result is that the Galois is not 
solvable. It is thus not abelian and the Hamiltonian ( 1161) is not integrable 
in the meromorphic function category on the chosen, six-dimensional phase 
space. To check the Liouville integrability in a concrete case like ( 152"]) one can 
also use an implementation of the Kovacic algorithm in symbolic manipulation 
programs, e.g. kovacicsols from Maple 12 which returns a list of Liouvillian 
solutions if they exist and the empty set in the opposite case. The occurrence 
of the latter case is thus a proof of the nonintegrability of the full system. 

6 Conclusions and outlook 

We proved in an analytical way the non-integrability of a specific quadratic 
Hamilton function defined on the Lie algebra SU3 obtained as a classical limit 
of a quantum Hamiltonian. The motivation was thoroughly presented in the 
Introduction and Section [3j here we want to conclude that the achieved result 
fills, at least partially, a gap in the reasoning of [8] . 

We would like also to highlight some novelties of our investigation. Since 
the birth of the Morales- Ramis theory there has been many successful at- 
tempts to apply it to concrete physical situations [231 EH ES1 ES]- The in- 
vestigated system were, usually, of the standard type with Hamilton function 
of the form of a sum of the kinetic and potential energies, the kinetic en- 
ergy being a quadratic form in the canonical momentum variables, defined 
in a topologically simple phase space. This is not the case for the Hamilton 
function treated in our paper. It is a quadratic polynomial in non-canonical 
variables on a compact symplectic manifold. 

The presented reasoning can be applied to other symplectic leaves, other 
Hamilton functions on 5113 (or other Lie algebras). The desired result would be 
a classification of such Hamilton functions with respect to their integrability. 
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